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Under  Contract  No.  F49620-C-89-0038,  NTNF/NORSAR  is  conducting 
research  within  a  wide  range  of  subjects  relevant  to  seismic  monitoring.  The 
emphasis  of  the  research  program  is  on  developing  and  assessing  methods 
for  processing  of  data  recorded  by  networks  of  small-aperture  arrays  and  3- 
component  stations,  for  events  both  at  regional  and  teleseismic  distances.  In 
addition,  more  general  seismological  research  topics  are  addressed. 

Each  quarterly  technical  report  under  this  contract  presents  one  or  several 
separate  investigations  addressing  specific  problems  within  the  scope  of  the 
statement  of  v.  ork.  Summaries  of  the  research  efforts  within  the  program  as  a 
whole  are  given  in  annual  technical  reports. 

This  Scientific  Report  No.  6  presents  a  manuscript  entitled  “Application 
of  Regional  Arrays  in  Seismic  Verification”,  by  S.  Mykkeltveit,  F.  Ringdal,  T. 
Kvaerna  and  R.W.  Alewine. 
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Application  of  Regional  Arrays  in  Seismic  Verification 

Research 

SvEIN  MYKKBLTVEIT*,  FrODE  RINGDAL*,  TORMOD  KViGRNA*  and" 
Ralph  W.  Alewine** 

*  NTNF/NORSAR,  P.O.  Box  51,  N-2007  Kjeller,  Norway 
*+  DARPA/NMRO,  1400  Wilson  Blvd.,  Arlington,  Virginia  22209,  USA 


Abstract 


The  paper  gives  an  account  of  the  work  related  to  the  development  of  the  NORESS 
concept  of  a  regional  array.  The  array  design  considerations  and  objectives  are  reviewed, 
and  a  description  is  given  of  the  NORESS  and  ARCESS  array  facilities  in  Norway  with 
their  field  installations,  data  transmission  lines  and  data  receiving  center  functions. 

The  automatic  signal  detection  processing  of  NORESS  data  applies  multiple  narrow- 
band  frequency  filters  in  parallel  and  forms  array  beams  from  selected  subgeometries.  The 
detection  algorithm  is  based  on  computing  the  STA/LTA  ratio  for  each  beam  individually, 
and  a  detection  is  declared  whenever  this  ratio  exceeds  a  given  threshold.  It  is  explained 
how  the  beam  deployment  and  the  individual  threshold  values  can  be  tuned  to  ensure 
that  the  interesting  phase  arrivals  are  not  missed,  but  at  the  same  time  avoiding  coda 
detections. 


For  each  detected  signal,  frequency-wavenumber  analysis  is  invoked  to  determine  ar¬ 
rival  azimuth  and  apparent  velocity.  Currently,  a  broad  band  estimator  is  used,  and  it 
is  demonstrated  that  use  of  this  algorithm  increases  the  stability  of  the  azimuth  and  ap¬ 
parent  velocity  estimates,  relative  to  narrow  band  methods.  Local  and  regional  events 


are  automatically  located  on  the  basis  of  identification  and  asrociation  of  P-  and  S-\ 


arrivals.  The  uncertainty  in  the  arrival  azimuth  is  the  limiting  factor  in  accurately  de¬ 
termining  single-array  event  locations,  and  it  is  shown  that  this  uncertainty  is  as  large  as 
10°-15°  for  Pn  phases  from  certain  regions. 
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In  order  to  father  investigate  the  potential  of  the  NORKSS  concept,  work  was  initiated 
towards  installing  a  network  of  regional  arrays  in  the  northern  Europe  area.  This  involved 
the  deployment  of  the  ARCESS  array  in  northern  Norway,  and  the  installation  of  the 
FINESA  array  in  Finland  in  cooperation  with  the  University  of  Helsinki.  Data  from  these 
three  arrays  have  been  used  jointly  in  a  location  estimation  scheme.  It  is  shown  that  for 
events  in  the  Fennoscandian  region  of  magnitude  typically  around  2.5  and  for  which  at 
least  one  phase  is  detected  by  each  array,  location  estimates  can  be  obtained  automatically 
that  deviate  from  published  network  locations  by  only  16  km  on  the  average. 

In  the  future,  it  is  anticipated  that  additional  arrays  and  single  stations  in  the  northern 
Europe  area  will  contribute  real-time  data  to  NORSAR  for  analysis  jointly  with  existing 
arrays.  The  first  additional  data  to  become  available  will  be  from  the  GERESS  array, 
which  will  be  established  in  the  Federal  Republic  of  Germany  in  1990.  Future  perspectives 
also  include  the  use  of  expert  system  technology  in  the  data  analysis,  and  the  IMS  system 
already  in  operation  represents  the  initial  attempt  in  this  regard.  A  summary  is  given  of 
problem  areas  where  further  work  is  needed  in  order  to  fully  exploit  the  regional  array 
concept. 


Introduction 

The  suggestion  to  use  seismic  arrays  in  order  to  detect,  locate  and  identify  low- 
magnitude  events  for  the  purpose  of  verifying  compliance  with  nuclear  testing  treaties 
dates  back  to  the  Geneva  Conference  of  experts  in  1958.  In  the  1960s  and  70s,  seismic 
arrays  were  established  in  several  countries  around  the  world,  and  the  arrays  were  gener¬ 
ally  designed  for  optimum  detection  capabilities  for  events  at  teleseismic  distances.  The 
most  ambitious  undertaking  in  this  regard  was  the  deployment  of  the  LAS  A  array  in 
the  United  States  (525  short  period  seismometers  over  an  aperture  of  200  km)  and  the 
NORSAR  array  in  Norway  (132  short  period  seismometers,  array  aperture  100  km).  Over 
the  years,  operation  of  teleseismic  arrays  has  testified  to  their  excellent  performance  in 
detecting  weak  arrivals,  as  well  as  their  ability  to  estimate  the  direction  and  apparent 
velocity  of  incoming  signals.  A  detailed  review  of  these  developments  is  given  by  Ringdal 
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and  Ilusebye  (1982). 

The  trilateral  negotiations  during  1977-80  on  a  Comprehensive  Test  Ban  Treaty  prompt¬ 
ed  a  shift  of  interest  from  observations  made  at  telcscismic  distances  to  wave  propagation 
in  the  regional  regime  (up  to  2-3000  km).  It  is  in  this  context  that  experiments  were 
initiated  in  Norway  in  1979  towards  the  development  of  a  “prototype”  regional  array, 
suitable  for  monitoring  of  low-level  seismic  activity  within  regional  distance  range.  It  was 
anticipated  that  this  work  would  be  aided  and  facilitated  by  the  experience  and  knowledge 
gained  from  10  years  of  operation  of  the  NORSAR  array,  but  it  was  also  realized  that  new 
experimental  data  had  to  be  obtained  to  adequately  design  an  array  with  the  desirable 
performance  for  regional  seismic  phases. 

The  purpose  of  this  paper  is  both  to  offer  an  overview  of  the  work  conducted  since  1979 
relating  to  the  development  of  the  NORESS  and  ARCESS  arrays,  and  at  the  same  time  to 
give  an  assessment  of  the  capbilities  of  these  arrays.  The  regional  array  program  with  its 
associated  research  activities  from  its  inception  in  1979  has  grown  to  become  significant 
in  both  size  and  diversity.  A  complete  review  of  these  developments  is  beyond  the  scope 
of  this  paper.  We  will  instead  focus  upon  those  aspects  of  these  developments  that  in 
our  judgement  are  the  most  important  ones  in  the  seismic  monitoring  context.  This  is 
done  by  firstly  reviewing  the  considerations  that  went  into  the  NORESS  design  efforts, 
before  details  are  given  on  the  array  installations  and  the  NORESS  and  ARCESS  field 
sites.  The  various  steps  in  the  automatic  data  processing  are  described,  and  the  individual 
performance  of  each  of  the  arrays  in  detecting  and  locating  events  at  regional  distances  is 
assessed.  Thereupon,  we  consider  the  capabilities  of  a  network  of  NORESS-type  arrays.  In 
making  our  assessments,  we  summarize  important  findings  available  in  the  literature,  and 
supplement  these  with  hitherto  unpublished  results  from  our  own  recent  research.  Finally, 
the  results  obtained  during  these  past  ten  years  are  discussed,  and  some  perspectives  for 
the  future  are  given. 


NORESS  Design  Considerations 

The  desirable  characteristics  of  a  prospective  prototype  regional  array  were  formulated 
at  the  outset  of  the  experiments  initiated  in  1979:  The  array  should  be  designed  for 
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optimum  detection  of  regional  seismic  signals,  and  it  should  provide  suflicient  resolution 
to  reliably  estimate  the  apparent  velocity  and  azimuth  of  all  such  signals.  Furthermore, 
it  was  dear  that  the  desirable  performance  of  the  array  with  respect  to  signal  detection 
and  characterization  would  need  to  be  obtained  over  the  wide  range  of  frequencies  typical 
of  regional  wave  propagation.  These  requirements  can  be  formulated  in  a  more  technical 
language,  as  follows: 

•  The  array  should  provide  close  to  optimum  gain  by  beamforming  for  the  phases  and 
frequencies  characteristic  of  regional  wave  propagation. 

•  The  array  geometry  should  be  symmetric  in  order  to  offer  equal  capabilities  for 
signals  from  all  directions,  and  the  response  pattern  should  have  a  narrow  main  lobe 
and  small  side  lobes. 

It  was  dear  at  the  outset  that  data  from  the  NORSAR  teleseismic  array  could  not 
be  used  to  infer  an  optimum  array  configuration  for  a  regional  array,  as  regional  signals 
recorded  at  NORSAR  with  its  minimum  sensor  separation  of  the  order  of  2.5  km  are 
spatially  aliased  and  furthermore  do  not  correlate  well  even  across  the  10  km  aperture 
subarrays.  The  main  emphasis  in  the  initial  experiments  was  therefore  placed  on  deriving 
signal  and  noise  correlation  curves  for  various  frequency  bands  for  intersensor  separations 
in  the  distance  interval  0-2  km.  Such  curves  can  be  used  to  express  the  beamforming  gain 
G  via  the  formula: 

N  N 

G"  =  E  Cij  /  Yi  pii  (1) 

«.i=i  «.;=i 

where  C,j  is  the  signal  correlation  between  sensors  i  and  j,  pij  is  the  corresponding  noise 
correlation,  and  N  is  the  number  of  sensors.  The  provisional  configuration  deployed  in 
1979  comprised  only  6  instruments  unevenly  spaced  within  an  aperture  of  2  km.  Still, 
signal  and  noise  correlation  curves  were  obtained  that  possessed  most  of  the  characteristic 
features  and  thus  qualitatively  resembled  the  curves  derived  later  on  from  configurations 
comprising  many  more  sensors.  Examples  of  such  curves,  derived  from  the  eventual  25- 
sensor  NORESS  geometry,  are  shown  in  Fig.  1. 

Analytical  representations  of  the  very  early  versions  of  the  signal  and  noise  correlation 
curves  were  used  by  Mykkeltveit  ei  al  (1983)  in  equation  (l)  to  find  geometries  that  max- 
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iinized  the  array  gain  G.  It  was  demonstrated  in  that  paper  that  optimized  geometries 
could  be  obtained  that  were  associated  with  theoretical  gains  well  in  excess  of  the  standard 
y/N  gain  by  utilizing  negative  minima  in  the  observed  noise  correlation  curves  (see  Fig. 
1).  Such  optimum  geometries,  however,  tended  to  be  rather  “peaked”  in  their  frequency 
response,  :.e.,  a  very  high  gain  at  one  particular  frequency  was  g.enerally  accompanied 
by  low  gains  at  other  frequencies.  The  optimized  geometries  were  characterized  by  one 
particular  intersensor  spacing  being  represented  as  many  times  as  possible  in  the  geome¬ 
try.  This  distance  reflected  the  separation  for  which  the  noise  correlation  curve  attained 
its  minimum,  for  a  given  frequency  interval.  For  optimization  explicitly  taking  several 
frequency  bands  into  consideration  (e.g.,  by  giving  equal  weight,  in  the  gain  expression,  to 
each  of  five  different  frequency  bands),  again  one  single  intermediate  frequency  dominated 
the  geometry.  At  this  stage,  it  was  realized  that  it  would  be  difficult  to  arrive  at  con¬ 
figurations  with  a  sufficiently  broad  frequency  response,  following  this  strategy.  Instead, 
the  approach  of  estimating  the  gain  via  equation  (1)  was  pursued  in  combination  with 
design  ideas  set  forth  by  Followill  and  Harris  (1983).  They  proposed  a  geometry  based  on 
concentric  rings  spaced  at  log-periodic  intervals  in  radius  iZ,  according  to  the  relation: 

R  =  limin' n  =  0,l,2,3  (2) 

The  geometry  of  the  NORESS  array  deployed  in  1984  is  a  realization  of  (2),  with  Rmin 
=  150  m  and  a  =  2.15.  Additional  details  on  how  the  partly  conflicting  demands  made 
on  array  performance  were  balanced  by  the  adoption  of  this  configuration  can  be  found  in 
Mykkeltveit  (1985). 

The  NORESS  array  configuration  is  shown  in  Fig.  2.  There  is  a  center  element  denoted 
AO,  3  elements  in  the  innermost  ring  (A-ring,  nominal  radius  150  m),  5  elements  in  the 
B-ring  (radius  323  m),  7  elements  in  the  C-ring  (radius  693  m),  and  9  elements  in  the  D- 
ring  (radius  1491  m,  giving  an  array  aperture  of  approximately  3  km).  The  short  period 
stations  at  AO,  C2,  C4  and  C7  are  equipped  with  3-component  instruments.  It  is  readily 
.seen  from  Fig.  2  that  theie  la  a  subbtaittia]  range  of  intcrsensor  separations  present  in  the 
NORESS  geometry,  and  that  it  offers  a  possibility  of  using  widely  different  subgeometries 
for  different  signal  frequencies. 


Description  of  NORESS  and  ARCESS 


The  NORESS  array  was  installed  in  southeastern  Norway  in  the  fall  of  1984  (See  Fig. 
13  for  location;  array  center  coordinates  are  60.735“N,  11.541°E)  as  a  joint  undertaking 
between  NORSAR,  the  Defense  Advanced  Research  Projects  Agency  and  the  Sandia  Na¬ 
tional  Laboratories.  The  array  site  is  in  a  wooded  area  with  relatively  low  population 
density.  There  is  no  human  activity  within  the  array  except  occasional  ferest  work  by 
the  landowner.  Competent  bedrock  with  P  wave  velocities  of  the  order  of  5.S-6.0  km/s 
is  found  either  at  the  surface,  or  underneath  a  layer  of  soil  of  tliickness  up  to  a  few  me¬ 
ters.  The  rock  is  of  Precambrian  age  and  is  composed  of  gneisses  and  gabbro.  A  seismic 
reflection  profile  running  north-south  slightly  east  of  the  array  center  showed  strong  in¬ 
dications  of  a  dipping  reflector  intersecting  the  surface  in  the  southeastern  part  of  the 
array  (Mykkeltveit,  1987).  Tomographic  mapping  of  the  valocity  structure  of  the  upper 
few  kilometers  beneath  NORESS  has  been  attempted  by  Ruud  and  Husebye  (1990). 

All  NORESS  short  period  instruments  are  placed  in  shallow  vaults  on  concrete  pads 
anchored  to  the  bedrock.  The  seismometer  housing  is  a  fiberglass  construction  sealed  to 
the  concrete  pad  to  prevent  water  leakage.  The  3-component  broad  band  seismometer 
at  the  array  center  is  deployed  in  a  60  m  deep  borehole.  All  data  from  the  vaults  and 
borehole  are  transmitted  to  the  hub  building  at  the  array  center  via  trenched  fiber  optic 
cables,  which  are  used  since  they  are  immune  to  electrical  disturbances  from,  e.ff.,  nearby 
power  lines  and  lightning.  Each  seismometer  site  is  powered  via  buried  cables  from  the 
hub  building. 

The  NORESS  short  period  instruments  are  of  type  GS-13,  and  the  broad  band  borehole 
seismometer  is  a  KS-36000  instrument.  The  conversion  of  data  from  analog  to  digital  form 
takes  place  in  the  vaults  and  in  the  borehole.  A  16  bits  A/D  converter  is  used,  with  2  of  the 
16  bits  used  for  gain  ranging  in  steps  of  1,  8,  32  and  128.  Short  period  data  are  digitized 
at  a  rate  of  40  Hz,  whereas  data  from  the  broad  band  instrument  are  sampled  both  at  1  Hz 
for  a  long  period  band  and  at  10  Hz  for  an  intermediate  period  band.  The  system  transfer 
functions  for  the  various  passbands  are  shown  in  Fig.  3.  The  high  frequency  station  that 
was  integrated  in  NORESS  in  1985  uses  the  analog  output  from  the  3-component  short 
period  instrument  at  site  AO,  and  digitizes  this  data  at  a  rate  of  125  Hz,  using  a  24  bits 
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A/D  converter.  A  detailed  description  of  this  system  is  given  by  lUngdal  ct  al  (1990). 

Tl\c  AR.CESS  array  installed  in  northern  Norway  in  the  fall  of  1987  (array  center 
coordinates  are  G9.53'1°N,  25.511°E;  see  location  in  Fig.  13)  represented  the,  first  step 
towards  a  network  of  NORESS-type  arrays.  ARCESS  is  located  at  a  distance  of  1174  km 
from  NORESS  in  an  area  of  low  population  density  and  little  or  no  industrial  activity. 
The  array  sensors  are  deployed  on  gabbro,  which  is  mostly  exposed  since  the  soil  cover  is 
nonexistent  or  very  thin  (up  to  0.5  m).  The  short  period  seismometers  are  located  in  drums 
placed  on  the  surface  and  covered  with  turf  and  moss,  and  the  broad  band  instrument  is 
in  a  50  m  deep  borehole.  Otherwise,  the  ARCESS  field  installation  closely  resembles  that 
of  NORESS:  The  geometries  are  nearly  the  same  (deviations  of  the  order  of  a  few  tens 
of  meters  in  relative  sensor  positions  exist  due  to  adjustment  to  local  terrain),  and  the 
seismic  system  with  sensors  and  other  electronic  components  are  identical. 

All  data  from  the  field  installations  of  NORESS  and  ARCESS  are  collected  by  the 
hub  processors  at  the  central  sites  and  transmitted  in  real  time  to  the  NORSAR  data 
processing  center  at  KJeller.  NORESS  data  are  transmitted  over  a  64  Kbits/s  land  line, 
whereas  a  domestic  satellite  link  with  the  same  capacity  is  used  to  transmit  ARCESS 
data.  At  Kjeller,  the  data  are  acquired  on  cyclic  disk  buffers  that  hold  72  hours  of  data 
for  each  array,  processed,  and  permanently  archived  on  magnetic  tapes  (on  8  mm  video 
cassettes  from  February  1990).  In  this  way,  data  from  the  last  3  days  can  be  directly 
accessed  from  disk,  whereas  any  data  can  be  retrieved  from  the  archive.  NORESS  and 
ARCESS  data  contain  a  substantial  amount  of  environmental  information  (temperatures, 
humidities,  wind  speed  and  direction),  as  well  as  state-of-health  and  instrument  calibration 
data,  which  are  being  analyzed  automatically  at  Kjeller  to  assist  in  detecting  system 
malfunction.  The  total  amount  of  data  generated  by  each  array  per  day  and  s*^orcd  in  the 
archives  is  approximately  400  Mbytes. 

Automatic  Data  Processing 

With  the  amount  of  data  received  continuously  from  arrays  like  NORESS,  it  is  of 
paramount  importance  that  reliable  schemes  for  fully  automatic  event  detection  and  lo¬ 
cation  in  near  real  time  be  developed  and  implemented.  The  development  of  such  algo- 


rithms  and  procedures  for  NORESS  went  in  parallel  with  the  array  design  work  in  the 
early  1980s.  The  RONAPP  (Regional  ON-line  Array  Processing  Package)  code,  described 
in  Mykkcltveit  and  Bungum  (1984)  was  the  result  of  these  efforts.  In  short,  RONAPP 
detects  phase  arrivals  using  an  STA/LTA-detector  applied  to  a  number  of  beams,  esti¬ 
mates  arrival  azimuth  and  apparent  velocity  for  the  signals  detected,  and  associates  P- 
and  S-arrivals  for  event  location.  Fig.  4  serves  as  an  illustration  of  this  procedure.  It 
shows  NORESS  data  for  a  presumed  nuclear  explosion  in  the  White  Sea  region  of  the 
USSR  on  July  18,  1985,  at  a  distance  of  1550  km  from  NORESS.  RONAPP  detected  Pn, 
Sn  and  Lg  arrivals  from  this  event,  as  indicated  in  the  figure.  The  information  on  arrival 
azimuth  and  apparent  velocity  derived  by  computing  frequency-wavenumber  spectra  for 
short  data  segments  around  the  arrival  times  is  used  together  with  standard  travel  time 
tables  to  identify  the  phases  and  locate  the  event. 

The  RONAPP  processing  package  has  evolved  through  several  generations  into  the 
version  currently  in  use.  The  main  underlying  ideas,  however,  remain  the  same  as  those 
described  in  detail  in  Mykkeltveit  and  Bungum  (1984).  In  the  following  description  of  the 
current  version  of  the  package  we  therefore  focus  on  those  aspects  that  have  undergone 
changes  since  1984  in  light  of  experience  gained. 

Signal  detection 

The  NORESS  and  ARCESS  detection  processing  is  similar  to  that  of  most  other  arrays: 
to  enhance  weak  signals,  a  number  of  filtered  beams  (steered  towards  various  hypothetical 
epicentral  locations)  are  computed  in  real  time  and  subjected  to  a  conventional  STA/LTA 
detector.  Whenever  the  STA/LTA  ratio  exceeds  a  preset  threshold,  a  detection  is  declared. 
When  the  threshold  is  exceeded  for  several  beams  within  a  4  s  long  window,  only  one 
detection  is  declared  and  it  is  attributed  to  the  beam  with  the  highest  STA/LTA  ratio. 

The  real  challenge  in  the  design  of  the  detector  is  the  selection  of  a  proper  beam 
deployment  with  associated  detection  thresholds.  Fig.  5  serves  to  illustrate  some  of  the 
considerations  involved.  It  shows  NORESS  data  from  a  small  event  in  western  Norway, 
at  an  epicentral  distance  of  approximately  350  km.  The  regional  phases  Pn,  Pg,  Sn  and 
Lg  are  clearly  seen  after  enhancement  by  filtering  the  data  in  widely  different  frequency 
bands  and  also  displaying  different  components  of  the  3-axis  station  at  site  AO.  It  is  found 
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for  this  particular  event  and  many  other  events  as  well  that  the  Pn  and  Pg  phases  are 
best  observed  on  vertical  channels,  but  in  different  frequency  bands  (here:  10-16  Hz  for 
Pn  and  3.5--5.5  Hz  for  Pg).  'J'he  Lg  phase  stands  out  clearly  in  the  1-2  Hz  band,  also 
on  vertical  channels.  The  onset  of  the  Sn  phase,  however,  is  very  often  found  to  have  an 
impulsive  character  on  the  horizontal  channels,  and  in  a  relatively  high  filter  band  (here: 
5-8  Hz).  This  examples  shows  that  a  beam  deployment  designed  for  detection  of  regional 
phases  must  Include  several  beams  filtered  within  different  narrow  frequency  bands,  for 
each  steering  direction.  It  is  also  important  to  utilize  the  horizontal  components  of  the 
3-component  stations  for  detection  of  phases  like  Sn. 

The  NORESS  beam  deployment  in  use  since  13  April  1989  is  given  in  Table  1.  It  is 
composed  of  76  beams,  out  of  which  66  are  conventional,  coherent  ones.  These  66  beams 
are  aimed  at  detection  of  P-phases  at  all  frequencies  and  are  designed  and  deployed  in 
accordance  with  the  criterion  that  the  gain  loss  due  to  missteering  should  be  less  than 
3  dB  for  any  signal  from  any  direction,  arriving  at  NORESS  with  an  apparent  velocity 
above  6.0  km/s.  The  subconfiguration  defined  for  each  of  these  beams  is  the  one  that  has 
been  found  to  provide  the  best  SNR  gain  for  the  frequency  band  in  question,  see  Kvaerna 
(1989).  Incoherent  beams  are  particularly  suited  for  detection  of  secondary  phases,  which 
are  often  of  an  emergent  nature.  Ten  such  beams  are  included  in  the  NORESS  beam  set, 
and  are  specifically  aimed  at  detecting  Sn  (beams  NIlOl-04)  and  Lg  (beams  NVOl-06) 
arrivals.  Details  on  incoherent  beamforming  are  given  in  Ringdal  et  al  (1975). 

The  determination  of  detection  thresholds  needs  special  attention,  and  one  important 
aspect  is  the  balancing  of  thresholds  between  the  coherent  and  incoherent  beams.  The 
thresholds  given  in  Table  1  were  determined  on  the  basis  of  operational  experience  as  well 
as  theoretical  considerations  on  false  alarm  rates.  Details  of  how  this  was  done  are  given 
in  Kvaerna  et  al  (1987). 

The  beam  set  used  at  NORESS  during  1  January  1985-13  April  1989  comprised  20 
beams  (see  Ringdal,  1990),  out  of  which  3  were  incoherent,  and  with  individual  thresholds 
very  similar  to  the  ones  of  the  new  beam  set.  The  ‘old’  beam  deployment  resulted  in 
approximately  50,000  detections  per  year,  and  the  new  and  more  extensive  beam  set  has 
caused  an  increase  in  the  number  of  detections,  mainly  due  to  an  improved  beam  coverage 
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for  high  frequencies.  The  ARCESS  beam  set  is,  since  1989,  identical  to  that  of  NORESS, 
except  for  the  steering  parameters  of  the  special  beams  directed  towards  Novaya  Zemlya 
and  Semipalatinsk.  The  number  of  detections  on  ARCESS  exceeds  that  of  NORESS;  see 
Bratt  et  al  (1990). 

Estimation  of  signal  attributes 

Following  the  detection  of  a  signal,  best  estimates  of  the  arrival  time,  signal  frequency 
and  amplitude  are  obtained  as  outlined  in  Mykkeltveit  and  Bungum  (1984).  A  frequency- 
wavenumber  (f-k)  spectrim  is  then  computed  for  a  3  s  long  time  window,  starting  0.5  s 
before  the  detection  time.  Initially,  the  narrow-band  f-k  method  described  in  Capon  (1969) 
was  used,  but  in  1989  a  wide-band  method  described  by  Kvaerna  and  Doornbos  (1986)  was 
adopted.  This  analysis  gives  estimates  of  the  arrival  aizimuth  and  apparent  velocity  of  the 
detected  signal.  The  signal  is  classified  as  P  or  S,  corresponding  to  the  apparent  velocity 
being  above  or  below  6.0  km/s,  respectively.  In  addition,  various  attributes  characterizing 
the  particle  motions  are  extracted  from  the  3-component  stations,  in  accordance  with  a 
scheme  devised  by  Jurkevics  (1988).  These  attributes  can  be  useful  to  distinguish  Pn 
phases  from  Pg,  and  Sn  phases  from  Lg,  although  consistent  separation  is  difficult  to 
achieve. 

Phase  association  and  event  location 

Simple  rules  and  procedures  are  used  to  associate  regional  phases  and  determine  event 
locations,  based  on  the  signal  attributes  estimated  during  the  post-detection  processing. 
These  rules  and  procedures  are  given  explicitly  in  Mykkeltveit  and  Bungum  (1984).  Mod¬ 
ifications  since  then  have  e.ssentially  amounted  to  the  following: 

•  Polarization  information  from  the  3-coinponent  stations  is  used  in  an  attempt  to 
determine  whether  a  P  arrival  is  Pn  or  Pg,  and  whether  an  S-arrival  is  Sn  or  Lg. 

•  All  regional  signals  with  a  phase  assignment  are  used  in  the  event  location  (previ¬ 
ously,  only  Pn  and  Lg  were  used).  Events  aie  located  using  the  TTAZLOC-algorlthm 
by  Bratt  and  Bache  (1988). 

•  Ph.ises  are  associated  if  their  arrival  azimuths  deviate  by  less  than  30°  (previously 
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20°),  and  arrivals  with  apparent  velocities  exceeding  12  km/s  are  not  considered 
regional  and  are  not  used  in  event  location. 


Experience  from  more  than  5  years  of  continuous  operation  of  NORESS  (and  more 
than  2  years  with  .A.RCESS)  has  shown  that  the  procedure  described  in  this  chapter  is 
successful  in  completely  automatically  detecting  and  locating  with  a  reasonable  accuracy  a 
substantial  number  (typically  10-20  per  day,  depending  on  day-of-week)  of  small  regional 
events.  In  particular,  we  want  to  emphasize  that  it  has  proved  to  be  possible  to  design 
the  detector,  with  its  beam  definitions,  associated  thresholds,  etc.,  in  such  a  way  that 
detections  are  declared  in  the  routine  processing  for  almost  all  regional  phases  that  can  be 
caught  by  the  human  eye  (including  those  that  can  only  be  seen  on  the  beams),  without 
causing  a  high  number  of  irrelevant  detections,  e.g.,  in  the  codas  following  the  phase 
onsets.  Too  many  detections  will  generally  cause  degradation  of  the  performance  of  the 
phase  association  and  event  location  steps.  Fig.  6  serves  to  illustrate  this  point.  It  shows 
the  final  output  from  the  automatic  processing  of  NORESS  data  from  a  mining  explosion 
in  Estonia,  at  an  epicentral  distance  of  795  km.  The  Pn,  Sn  and  Lg  phases  have  been 
detected  and  identified  as  such  and  subsequently  used  in  the  event  location  step.  The 
Pn  arrival  was  detected  on  the  N073  beam  with  an  STA/LTA  ratio  of  5.3.  This  beam  is 
shown  as  trace  no.  2  from  the  bottom  of  the  plot.  It  has  a  steering  azimuth  of  90°,  which 
is  close  to  the  azimuth  of  96.5°  estimated  by  the  f-k  analysis  for  this  phase.  The  Sn  phase 
was  detected  on  beam  NH02  (this  beam  is  not  shown),  and  the  Lg  phase  on  beam  NV04, 
which  is  displayed  in  the  form  of  a  coherent  beam  as  the  bottom  trace  in  Fig.  6.  Two 
additional  detections  in  the  Lg  coda  (indicated  by  small  arrows  above  the  upper  trace) 
show  Lg  type  phase  velocities  and  cause  no  difficulties  in  the  automatic  phase  association 
and  event  location  steps. 


In  addition  to  the  regional  arrivals,  many  teleseismic  signals  are  detected  by  NORESS 
and  ARCESS.  In  fact,  the  beam  deployment  in  Table  1  was  designed  also  with  the  tcleseis- 
mic  detection  performance  in  mind,  and  several  studies  (c.y.,  PJngdal,  1990)  have  testified 
to  the  excellent  capabilities  of  NORESS  and  ARCESS  in  this  regard.  Fluctuations  in  the 
seismic  noise  field  also  give  rise  to  many  detections  at  both  NORESS  and  ARCESS,  and 
these  are  usually  characterized  by  very  low  (Rayleigh  type)  apparent  velocities.  Kvaerna 
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(1990)  has  demonstrated  that  one  class  of  such  detections  at  NORESS  can  be  correlated 
with  the  waterflow  in  a  nearby  major  river. 

Capabilities  of  NORESS  and  ARCESS 

In  this  chapter,  we  will  review  and  assess  the  individual  capabilities  of  NORESS  and 
ARCESS  to  detect  and  characterize  regional  signals,  as^well  as  their  ability  to  locate 
regional  events. 

Signal  detection  capabilities 

Numerous  investigations  have  testified  to  the  excellent  capabilities  of  NORESS  and 
ARCESS  to  suppress  the  noise  and  thus  obtain  considerable  beamforming  gain.  NORESS 
noise  suppression  spectra  taken  hourly  during  a  one- week  period  in  July  1986  are  shown 
in  Fig.  7.  A  noise  suppression  spectrum  is  estimated  as  the  ratio  of  the  beam  power 
spectrum  to  the  average  power  spectrum,  taken  over  all  contributing  sensors.  For  the 
specific  subconfiguration  (AO,  C-  and  D-ring  sensors)  used  in  Fig.  7  the  noise  suppression 
is  particularly  effective  in  the  1.3-2.7  Hz  band,  where  the  suppression  exceeds  the  VW  level 
expected  for  uncorrelated  noise  by  up  to  6  dB.  The  trend  of  the  noise  suppression  curves  in 
this  frequency  interval  is  the  frequency-domain  manifestation  of  the  negative  correlations 
observed  in  the  curves  in  Fig.  1.  Additional  details  on  the  NORESS  noise  suppression 
capabilities  are  found  in  Mykkeltveit  et  al  (1990),  where  it  is  also  demonstrated  that 
these  capabilities  are  stable  over  time,  and  furthermore  show  no  strong  dependency  on 
the  actual  noise  level.  Noise  level  variations  at  NORESS  have  been  extensively  studied  by 
Fyen  (1990). 

Taking  also  signal  correlations  into  account,  Kvserna  (1989)  has  conducted  a  study 
to  determine  the  achievable  P  phase  SNR  gains  at  NORESS  for  various  frequencies.  He 
concludes  that  by  carefully  choosing  the  proper  subgeometry  at  different  frequencies,  gains 
exceeding  10  dB  can  be  consistently  achieved  over  almost  the  entire  band  0.5  to  10  Hz, 
with  the  highest  average  gains  of  12-14  dB  being  obtained  in  the  1-4  Hz  band. 

The  capability  of  NORESS  and  ARCESS  to  detect  regional  signals  can  be  inferred 
from  comparison  with  network  bulletins  for  the  region  under  study.  Adopting  tlie  regional 
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bulletin  published  by  the  University  of  Helsinki  as  a  reference,  we  have  associated  P- 
phases  detected  at  ARCESS  during  January  -  March  1988  with  events  reported  in  this 
bulletin  from  the  region  of  southern  Finland  and  the  surrounding  areas  of  the  western 
USSR,  with  epicentral  distances  relative  to  ARCESS  in  the  range  800-1200  km.  Fig.  8 
gives  a  histogram  of  the  number  of  reference  events  at  each  magnitude  with  the  events 
detected  by  ARCESS  marked  specially.  The  figure  further  contains  a  detection  probability 
curve  with  associated  confidence  limits,  estimated  by  the  maximum  likelihood  method  of 
Ringdal  (1975).  We  note  that  the  90  per  cent  P-wavc  detection  capability  in  the  region 
studied  is  close  to  Mf,  =  2.5.  This  can  be  compared  to  the  NORESS  threshold  of  Mi  = 
2.7  found  by  Ringdal  (1986)  for  detection  of  P  waves  for  events  in  the  same  region.  If  the 
criterion  is  relaxed  to  detection  of  either  a  P  or  an  S  arrival,  the  90  per  cent  detection 
threshold  is  approximately  0.2  magnitude  units  lower. 


The  Novaya  Zemlya  test  site  is  located  at  regional  distance  relative  to  both  ARCESS 
and  NORESS,  and  it  is  therefore  of  special  interest  to  take  a  closer  look  at  signals  from 
this  site.  Fig.  9  shows  data  for  an  m6  5.9  (NEIC)  explosion  on  December  4, 1988,  for  both 
ARCESS  and  NORESS.  For  both  arrays,  raw  (upper  trace)  and  filtered  data  (lower  trace) 
in  the  5-10  Hz  band  are  shown  for  the  sensor  at  Cl.  The  P  wave  was  detected  at  ARCESS 


with  an  STA/LTA  of  9136  (filter  band  3-5  Hz  for  the  detecting  beam),  and  at  NORESS 
with  an  STA/LTA  of  796  (filter  band  2.5-4.5  for  the  detecting  beam).  While  the  ARCESS 
detection  capability  is  definitely  better  than  that  of  NORESS  for  Novaya  Zemlya,  the 
difference  is  not  as  large  as  indicated  by  the  STA/LTA  values.  Thus,  the  reason  for  the 
low  value  for  NORESS  is  partly  that  the  first  P  arrival  has  a  low  amplitude  relative  to 
the  main  phase  that  arrives  about  4  seconds  later.  The  maximum  P  amplitude  is  in  fact 
similar  for  NORESS  and  ARCESS,  even  though  the  distance  to  NORESS  is  twice  as  large 
as  the  distance  to  ARCESS  (20°  and  10°,  respectively).  A  comparison  of  the  amplitudes  of 
the  data  filtered  in  the  5-10  Hz  band  shows  that  the  higher  frequencies  attenuate  rapidly 
with  distance.  The  secondary  phase  seen  at  ARCESS  is  Sn;  the  Lg  phase  is  absent,  most 


likely  due  to  structural  in  homogeneities  and  lateral  v3,riation  of  sediment  thicknesses  in 


the  Barents  Sea  (Baumgardt,  1990).  At  NORESS,  no  secondary  phases  are  readily  visible 


in  Fig.  9,  but  at  closer  examination  a  weak  Sn  phase  is  seen  in  a  low  frequency  passband. 
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Comprehensive  studies  on  spectral  characteristics  of  regional  phases  in  Fennoscandia 
have  been  carried  out  by  several  workers,  see,  e.g.,  Suteau-Henson  and  Baclie  (1988)  and 
Ringdal  el  al  (1990). 

Signal  esiimation  cajxibilities 

As  mentioned  in  the  previous  chapter,  a  broad  band  f-k  estimator  (Kvaerna  and  Doorn- 
bos,  1986)  has  replaced  the  narrow  band  f-k  analysis  technique  in  the  signal  estimation 
stage  of  the  online  data  processing  of  NORESS  and  ARCESS  data.  Before  this  was  done, 
the  performance  of  the  two  techniques  was  compared  on  several  different  data  sets.  One 
such  set  was  a  suite  of  10  chemical  explosions  at  a  dam  construction  site  (Blasj0)  in  south¬ 
ern  Norway  at  a  distance  of  301  km  from  NORESS.  The  two  analysis  techniques  were 
applied  to  the  Pn,  Sn  and  Lg  phases  from  the  10  events,  and  the  results  in  terms  of  esti¬ 
mated  arrival  azimuth  and  apparent  velocity  are  shown  in  Fig.  10.  It  is  readily  seen  that 
the  broad  band  method  offers  more  stable  results  than  the  narrow  band  technique  for  all 
phases.  The  separation  between  the  Sn  and  Lg  phases  is  quite  clear  for  this  data  set.  It 
is  difficult,  however,  in  general  to  separate  Sn  and  Lg  based  on  apparent  velocities  alone. 

Our  experience  with  this  broad  band  estimator  is  in  general  very  good,  and  it  appears 
to  be  the  best  tool  currently  available  for  estimation  of  arrival  azimuth  and  apparent 
velocity  from  NORESS  and  ARCESS  data.  Occasiondly  though,  it  is  observed  for  events 
with  known  epicenters  that  the  estimates  deviate  substantially  from  the  true  values,  and 
deviations  as  large  as  10°-15'’  are  observed  for  certain  source  regions  even  at  high  signal- 
to-noise  ratios.  These  arrivals  are  coming  in  off  azimuth  due  to  structural  inhomogeneities 
along  the  propagation  path,  and  the  estimates  just  reflect  true  propagation  characteristics. 

The  3-componenl  stations  at  NORESS  and  ARCESS  can  also  be  used  to  infer  slow¬ 
nesses.  These  are,  however,  less  accurate  than  those  obtained  by  the  f-k  analysis  using 
the  vertical  sensors  of  the  array,  particularly  at  lov/  values  of  SNR;  see  Harris  (1990)  and 
Henson  (1990).  In  addition,  3-componcnt  estimates  are  susceptible  to  surface  topogra¬ 
phy  (Odegaard  ci  cl,  1990).  However,  polarization  information  derived  from  3-component 
stations  can  aid  in  the  classification  of  the  phases  (Jepsen  and  Kennett,  1990). 

Event  location  capabilities 
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In  order  to  accurately  locate  regional  events  using  one-array  data,  both  the  epicentral 
distance  and  event  azimuth  must  be  reliably  determined. 

The  epicentral  distance  is  determined  on  the  basis  of  the  travel  time  dilTerence  between 
the  various  regional  phases.  Our  operational  experience  indicates  that  the  key  issue  in 
this  regard  is  the  proper  assignment  of  phase  types  (Sn,  Lg,  Rg)  for  the  S  waves  detected. 
When  only  one  S  wave  has  been  detected,  it  may  be  very  diflicult  to  decide  whether  it  is 
an  Sn  or  Lg  arrival.  An  illustration  of  this  is  provided  in  Fig.  11,  which  shows  the  relative 
strength  of  Sn  and  Lg  for  six  events  recorded  on  a  short  period  sensor  of  the  NORSAR 
array  (the  NORSAR  array  is  co-located  with  NORESS).  Fig.  11  shows  that  Lg  is  strongly 
attenuated  for  certain  propagation  paths  (as  also  observed  in  Fig.  9  for  the  Novaya  Zemlya 
to  ARCESS  and  NORESS  paths),  and  that  Sn  likewise  is  sometimes  not  seen  at  all  (event 
1).  It  is  well  known  that  these  differences  are  related  to  structural  inhomogeneities  along 
the  propagation  paths. 

A  misidentification  of  Sn  as  Lg  (or  vice-versa)  will  easily  cause  an  error  in  the  epi¬ 
central  distance  estimate  of  several  hundred  kilometers.  However,  when  phase  types  are 
correctly  assigned,  the  epicentral  distance  can  usually  be  determined  to  within  a  few  tens 
of  kilometers.  As  a  rule  of  thumb,  an  error  of  1  s  in  the  travel  time  difference  itween 
an  S  and  a  P  phase  results  in  an  error  of  6  km  in  the  epicentral  distance.  The  accuracy 
in  the  travel  time  difference  largely  depends  on  how  well  the  arrivad  time  of  the  S  phase 
can  be  determined,  and  emergent  arrivals  represent  a  problem  in  this  regard.  When  the 
phase  onsets  can  be  as  reliably  determined  as  shown  in  the  example  in  Fig.  6,  however,  the 
travel  time  tables  used  beconic  the  limiting  factor  for  accurate  determination  of  epicentral 
distance. 

Fig.  12  shows  azimuth  residuals  for  Pn  phases  detected  at  NORESS  during  1985- 
1988.  The  data  were  derived  as  follows:  The  NORESS  detection  lists  were  searched  for 
Pn  phases  that  could  be  associated  with  regional  events  reported  in  the  Helsinki  bulletin 
for  this  4-yoar  period.  'I'hc  NORESS  detection  lists  provide  information  on  the  arrival 
azimuth  estimated  for  these  Pn  phases.  The  azimuth  residuals  (estimated  azimuth  minus 
‘true’  azimuth)  were  computed  and  averaged  on  a  grid  of  1°  x  2°  blocks  (north-south 
and  east-west,  respectively).  These  average  values  are  then  represented  as  .shown  in  the 
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figure.  For  some  of  the  blocks  in  the  Estonian- Leningrad  region  of  the  USbl^,  averaging  is 
done  over  several  hundred  azimuth  residuals.  Fig.  12  shows  that  tlie  azimuth  residuals  are 
moderate  (less  than  2.r>°  in  absolute  value)  for  large  geographical  arca.s,  but  there  arc  also 
areas  with  complicated  behavior,  like  the  Estonian-Leningrad  region.  We  note  that  the 
arrival  azimuths  derived  during  1985-88  were  computed  using  the  narrow  band  f-k  method, 
and  that  the  broad  band  method  currently  being  used  is  likely  to  cause  a  reduction  in  the 
absolute  values  of  the  residuals.  Anyway,  the  event  locations  will  inevitably  reflect  such 
arrival  azimuth  uncertainties,  and  some  kind  of  regional  calibration  will  be  needed,  e.g., 
by  taking  information  such  as  presented  in  Fig.  12  into  account  and  correcting  the  arrival 
azimuths  accordingly. 


Network  of  NORESS-type  Arrays 

We  have  seen  in  the  aforegoing  chapter  that  single-array  event  locations  may  deviate 
from  true  locations  by  typically  several  tens  of  kilometers  and  even  more,  if  phases  are 
incorrectly  assigned.  The  next  step  to  improve  the  event  location  capability  would  be 
to  install  a  network  of  NORESS-type  arrays.  This  effort  started  with  the  installation  of 
the  ARCESS  array  in  1987,  and  continued  with  the  deployment  of  the  somewhat  smaller 
FINESA  array  in  Finland  (Uski,  1990)  in  cooperation  with  the  University  of  Helsinki, 
and  the  GERESS  array  in  the  Federal  Republic  of  Germany  (Ilarjes,  1990).  Data  from 
FINESA  are  since  1  January  1990  available  in  real  time  at  the  NORSAR  data  processing 
center  at  Kjeller  and  are  being  processed  in  the  same  way  as  data  from  ARCESS  and 
NORESS.  (The  first  version  of  the  FINESA  array  was  installed  in  1985,  but  before  1 
January  1990  data  were  recorded  in  trigger  mode  and  on  tapes  at  the  array  site  only.) 
Data  from  GERESS  will  be  available  from  the  summer  of  1990,  and  will  also  be  transmitted 
to  Kjeller  for  online  processing. 

In  the  following,  we  report  on  the  findings  of  an  investigation  where  data  recorded  at 
FINESA  during  experimental  operation  in  March  1988  were  used  together  with  NORESS 
and  ARCESS  data  in  assessing  the  capabilities  of  this  three-array  network  in  locating 
events  in  the  Fennoscandian  region.  A  set  of  10  event,  for  which  there  was  at  least  one 
detected  phase  for  each  array,  was  selected  for  an  event  location  exiieriinent.  'J'hc  events 
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arc  listed  in  Table  l  and  shown  in  Fig.  13.  The  event  magnitudes  range  from  less  than  2.0 
to  3.2.  The  epiccntral  locations  for  the  10  events  are  taken  from  the  Helsinki  bulletin. 

The  continuous  processing  of  data  recorded  at  each  of  the  three  regional  a.rrays  provides 
(among  others)  estimates  of  arrival  times,  arrival  azimuths,  and  indication  of  phase  type, 
as  explained  in  the  previous  chapters.  These  parameters  together  with  the  associated 
uncertainties  were  used  as  i.iput  to  the  TTAZLOC  program  developed  by  Bratt  and  Bache 
(1988).  TTAZLOC  incorjroratcs  the  arrival  time  and  azimuth  data  into  a  generalizcd- 
invcr.se  location  estimation  .scheme,  and  can  be  a|>plied  to  both  .single  array  and  mnlti|)le- 
array  data.  Table  2  gives  the  results  of  the  location  experiment.  On  the  average,  the 
joint  three-array  locations  deviate  from  the  network  locations  published  in  the  Helsinki 
bulletin  by  16  km.  Two-array  and  one-array  locations  were  computed  for  all  combinations 
of  events  and  array  sub-networks,  also  using  the  TTAZLOC  aigorithm.  The  resulting 
average  deviations  from  the  network  solutions  are  26  and  68  km,  respectively.  Bratt  and 
Bache  (1988)  and  Bratt  (1990)  have  also  considered  event  mis'oeations  using  one-array 
and  two-array  data,  and  their  results  agree  well  with  those  reported  here. 

We  see  that  the  result  for  two-array  locations  represents  a  significant  improvement 
over  one-array  locations,  and  that  there  is  a  further  improvement  when  invoking  data 
from  three  arrays.  We  consider  the  results  reported  here  as  quite  promising,  when  taking 
the  following  into  account: 

•  The  arrival  times  used  were  those  determined  automatically  by  the  online  processing. 
It  is  conceivable  that  human  intervention  for  adjustment  of  arrival  times  and/or 
refinement  of  the  automatic  procedure  would  improve  the  location  estimates. 

•  Only  standard  travel  time  tables  for  the  phases  Pn,  Sn  and  Lg  were  used.  The 
introduction  of  regionalized  travel  time  tables  is  likely  to  result  in  improvements. 

•  Master  event  location  schemes  of  various  kinds  hold  considerable  promise  and  are 
expected  to  further  enhance  the  capabilities  of  accurately  locating  regional  events. 

Discussion  and  Pkkspbctivbs  koiitiib  Futurb 

This  paper  has  summarized  the  main  results  from  ten  years  of  research  at  NORSAR 
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on  regional  arrays  and  their  capabilities  in  detecting  and  characterizing  seismic  phases 
from  events  at  regional  distance.  Assessments  have  also  been  given  of  the  accuracy  in 
determination  of  event  locations  from  one-,  two-  and  three-array  data.  In  summary,  we 
have  found  that  P  waves  from  events  in  the  magnitude  range  2.5-2.7  at  distances  around 
1000  km  are  detected  with  a  probability  of  90  per  cent,  that  the  detector  can  be  tuned  so  as 
to  ensure  detection  also  of  emergent  S  wave  arrivals,  that  the  important  signal  attributes 
of  regional  phases  (like  arrival  azimuth,  apparent  velocity,  state  of  polarization,  etc.)  can 
be  adequately  characterized  using  a  broad  band  f-k  estimator  on  the  vertical  sensors  of  the 
array  in  combination  with  techniques  for  analysis  of  3-component  data,  that  the  regional 
phases  can  be  associated  using  a  set  of  simple  rules,  and  finally,  that  regional  events  can 
be  located  with  an  accuracy  of  the  order  of  15  km  when  data  from  a  network  of  3  arrays 
are  used.  It  is  noteworthy  that  these  results  are  all  obtained  from  fully  automatic  data 
processing,  t.e.,  without  human  intervention  of  any  kind. 

In  several  places  in  tliis  paper,  we  have  dealt  with  the  issue  of  formulating  rules  that 
can  be  used  in  the  automatic  data  processing.  Such  rules  are  generally  suitable  for  the 
application  of  AI  technologies  in  the  form  of  so-called  expert  systems.  While  several  such 
systems  have  been  designed  and  experimentally  applied  to  processing  of  regional  array 
data  (e.jf.,  Hiebert-Dodd,  1989;  Mason  et  al,  1988),  the  most  ambitious  development  so 
far  is  the  Intelligent  Monitoring  System  (IMS),  described  by  Bache  et  al  (1990).  The  aim 
of  the  IMS  is  to  automate  as  much  as  possible  the  seismic  data  interpretation  process, 
thus  taking  advantage  of  results  such  as  those  reported  in  this  paper.  It  is  currently 
being  developed  into  a  complete  system  for  integrated,  automatic  processing  of  data  from 
a  network  of  regional  arrays  and  single  stations.  The  first  version  of  the  IMS  has  been 
operated  at  the  Center  for  Seismic  Studies  in  Arlington,  Virginia,  and  at  NORSAR  since 
the  fall  of  1989.  Initial  results  for  the  performance  of  IMS  on  data  from  NORESS  and 
ARCESS  are  given  by  Bratt  et  al  (1990). 

The  deployment  of  .  Ivanced  regional  arrays,  and  the  associated  development  and 
implementation  of  automated  and  increasingly  powerful  data  processing  techniques  repre¬ 
sents  one  of  the  major  advances  in  the  field  of  seismic  monitoring  in  recent  years.  Arrays 
of  the  NORESS/ARCESS  design  have  demonstrated  a  capability  tc  lower  the  detection 
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threshold  by  more  than  0.5  magnitude  units  over  a  wide  range  of  signal  frequencies  rel¬ 
ative  to  traditional  seismic  stations.  Furthermore,  such  regional  arrays  provide  reliable 
phase  identification  and  azimuth  estunates  which  are  particularly  useful  in  locating  weak 
events  that  are  detected  by  only  a  few  stations.  In  a  seismic  monitoring  context,  these  are 
precisely  those  events  that  will  need  to  be  given  the  most  emphasis. 

Still,  the  process  of  exploiting  the  full  potential  of  regional  arrays  is  only  at  its  begin¬ 
ning.  Much  research  remains  to  be  done  in  developing  methods  for  integrated  processing 
of  data  from  a  regional  network  composed  of  arrays  and  3-componcnt  stations,  and  to 
assess  its  capabilities  in  a  seismic  monitoring  context. 

An  essential  aspect  will  be  the  further  development  of  advanced  processing  technol¬ 
ogy  to  handle  the  data  acquisition  and  quality  control,  data  base  organization,  automatic 
interpretation  and  interactive  computer  graphics  that  will  be  required  in  an  actual  moni¬ 
toring  situation.  From  a  seismological  point  of  view,  research  topics  to  be  addressed  will 
include: 

Signal  processing:  An  important  task  will  be  to  improve  current  methods  for  auto¬ 
matic  detection,  phase  identification  and  onset  time  determination,  with  particular  em¬ 
phasis  on  Sn  and  Lg  phases.  Methods  should  be  developed  to  exploit  the  potential  for 
array  processing  of  3-component  recordings  ofTered  by  NORESS-type  arrays,  and  to  im¬ 
prove  polarization  analysis  techniques  at  low  SNR.  There  is  a  need  to  address  the  problem 
of  how  to  automatically  resolve  multiple  events,  t.e.,  two  or  more  events  with  detected 
phases  that  are  intermixed  in  time  {e.g.,  multiple  mining  explosions). 

Regional  calibration:  An  important  aspect  is  the  development  of  region-specific 
corrections  for  travel-time  and  azimuth  anomalies.  Furthermore,  as  attenuation  charac¬ 
teristics  of  various  seismic  phase  types  (Pn,  Pg,  Sn,  Lg,  Rg)  are  highly  path-dependent, 
it  will  be  essential  to  accumulate  region-specific  knowledge  in  this  regard.  This  will  also 
contribute  to  the  establishment  of  consistent  regional  magnitude  scales,  which  should  be 
developed  and  calibrated  so  as  to  be  compatible  with  teleseismic  magnitude. 

Master  event  technique:  The  most  powerful  method  currently  available  to  obtain 
very  precise  epicenter  and  depth  estimates  is  joint  hypocentral  determination  using  well- 
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recorded  master  events.  This  concept  should  be  systematically  extended  to  more  general 
script-based  pattern  matching  techniques,  using  the  full  array  capabilities  and  incorporat¬ 
ing  both  ti’  le  domain  and  frequency  domain  features.  This  would  be  particularly  useful  to 
monitor  specific  mining  locations  and  test  sites,  and  should  be  accompanied  by  optimized 
beam  and  filter  settings  for  such  areas  of  special  interest.  The  method  for  continuous 
threshold  monitoring  (Ringdal  and  Kvaerna,  1989)  should  be  further  developed  and  ap¬ 
plied  in  such  a  context. 

Source  identification:  Although  efforts  to  develop  regional  discriminants  so  far  have 
met  with  little  success,  this  essential  research  should  be  pursued,  aiming  at  obtaining 
systematic  rules,  developed  on  a  region-specific  basis.  An  important  task  will  be  to  identify 
mining  shots,  and  the  very  promising  approach  of  characterizing  spectral  attributes  of 
ripple-fired  explosions  should  be  given  particular  attention.  The  use  of  regional  surface 
waves  (e.g.,  Rg)  in  source  identification  is  also  an  important  area  of  further  research. 

The  key  to  further  progress  in  this  field  would  appear  to  be  development  and  ap¬ 
plication  of  region-specific  knowledge,  both  for  the  purpose  of  detection,  location,  depth 
estimation  and  source  characterization.  Much  of  the  necessary  knowledge  can  only  be 
obtained  through  experimental  operation  over  an  extended  period  of  time.  The  regional 
array  network  now  being  developed  will  in  a  unique  way  contribute  to  the  establishment 
of  such  a  knowledge  base.  The  associated  data  management  and  processing  facilities  will 
ensure  the  availability  of  the  data  and  the  dissemination  of  the  accumulated  knowledge, 
and  will  thus  aid  in  future  efforts  directed  toward  these  research  goals. 
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Table  Captions 


Table  1.  NORESS  beam  deployment.  The  table  gives  name  of  beam,  steering  velocity 
(in  km/s),  steering  azimuth  (in  degrees),  filter  band  (in  IIz),  STA/LTA  threshold,  and 
subconfiguration  (in  terms  of  wliich  rings  are  included;  the  sensor  at  the  central  site  AO 
participates  in  all  beams).  The  NHOl-04  beams  are  incoherent  and  use  the  8  horizontal 
channels  (ne  in  the  table,  for  north-south  and  east-west)  of  the  stations  at  AO,  C2,  C4  and 
C7.  The  NVOl-06  beams  are  also  incoherent,  and  use  vertical  sensors  as  indicated.  The 
remaining  beams  are  conventional,  coherent  ones,  using  vertical  channels  only.  The  six 
coherent  beams  NH032-N037  are  identical  except  for  the  steering  azimuths,  which  have 
values  of  30,  90,  150,  210,  270  and  330  degrees,  respectively.  The  same  pattern  repeats 
for  other  coherent  beams  further  down  the  table,  and  is  indicated  by  an  asterix  in  the 
azimuth  column.  Four  special,  coherent  beams  are  steered  towards  each  of  the  test  sites 
at  Semipalatinsk  and  Novaya  Zemlya  (at  azimuths  of  80  and  30  degrees,  respectively). 

Table  2.  Results  from  TTAZLOC  location  experiments  using  data  from  NORESS, 
ARCESS  and  FINESA.  Epicentral  locations  are  given  as  reported  in  the  Helsinki  bul¬ 
letin  for  a  set  of  ten  regional  events  that  occurred  between  March  12  and  March  18  of 
1988.  The  table  gives  the  deviation  from  these  reference  locations,  as  inferred  from  the 
TTAZLOC  experiment  described  in  the  text. 
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Beam 

Velocity 

.'Vzimutli 

Filler  band 

Threshold 

Conf. 

Non 

99999.9 

0.0 

0.5-  1.5 

4.0 

D 

N021 

99999.9 

0.0 

1.0-  3.0 

4.0 

CD 

N031 

99999.9 

0.0 

1.5-  3.5 

4.0 

CD 

N032-037 

11.0 

1.5- 3.5 

4.0 

CD 

N038 

15.8 

80.0 

1.5- 3.5 

3.5 

CD 

N039 

10.0 

30.0 

1.5-  3.5 

3.5 

CD 

N041 

99999.9 

0.0 

2.0-  4.0 

4.0 

CD 

N042-047 

10.2 

♦ 

2.0-  4.0 

4.0 

CD 

N048 

15.8 

80.0 

2.0- 4.0 

3.5 

CD 

N049 

10.0 

30.0 

2.0- 4.0 

3.5 

CD 

N051 

99999.9 

0.0 

2.5-  4.5 

4.0 

BCD 

N052-057 

8.9 

♦ 

2.5-  4.5 

4.0 

BCD 

N058 

15.8 

80.0 

2.5-  4.5 

3.5 

BCD 

N059 

10.0 

30.0 

2.5-  4.5 

3.5 

BCD 

N061 

99999.9 

0.0 

3.0-  5.0 

4.0 

BCD 

N062-0G7 

10.5 

* 

3.0-  5.0 

4.0 

BCD 

N068 

15.8 

80.0 

3.0-  5.0 

3.5 

BCD 

N069 

10.0 

30.0 

3.0-  5.0 

3.5 

BCD 

N071 

99999.9 

0.0 

3.5-  5.5 

4.0 

BC 

N072-077 

11.1 

♦ 

3.5-  5.5 

4.0 

BC 

N081 

99999.9 

0.0 

o 

CO 

o 

4.0 

BC 

N082-087 

9.5 

♦ 

o 

CO 

b 

^.0 

BC 

N091 

99999.9 

0.0 

5.0-10.0 

4.5 

BC 

N092-097 

10.5 

♦ 

5.0-10.0 

4.5 

BC 

NlOl 

99999.9 

0.0 

8.0-16.0 

4.5 

AB 

N102-107 

9.9 

♦ 

8.0-16.0 

4.5 

AB 

NIIOI 

99999.9 

0.0 

2.0-  4.0 

2.4 

ne 

NH02 

99999.9 

0.0 

3.5-  5.5 

2.4 

ne 

NH03 

99999.9 

0.0 

5.0-10.0 

2.4 

ne 

NH04 

99999.9 

0.0 

8.0-16.0 

2.5 

ne 

NVOl 

99999.9 

0.0 

0.5-  1.5 

2.5 

D 

NV02 

99999.9 

0.0 

1.0-  2.0 

2.5 

C 

NV03 

99999.9 

0.0 

1.5-  2.5 

2.5 

C 

NV04 

99999.9 

0.0 

2.0-  3.0 

2.5 

C 

NV05 

99999.9 

0.0 

2.0-  4.0 

2.4 

C 

NV06 

99999.9 

0.0 

3.5-  5.5 

2.4 

n 

S-/ 

Table  1 
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Event 

No. 

Network 

Lat.  Long. 

Mag. 

Ml 

No.  of 
phases 
used 

3-airay 

‘error’ 

(km) 

Average 

2-array 

‘error’ 

(km) 

Average 

1 -array 
‘error’ 
(km) 

1 

67.1 

20.6 

<2 

5 

19 

31 

36 

2 

59.5 

26.5 

2.5 

8 

9 

8 

39 

3 

60.93 

29.19 

2.4 

6 

34 

34 

34 

4 

59.5 

25.0 

2.3 

8 

8 

23 

95 

5 

63.2 

27.8 

2.5 

6 

32 

31 

41 

6 

58.33 

10.93 

2.7 

7 

16 

24 

44 

7 

69.6 

29.9 

2.9 

8 

4 

13 

45 

8 

59.3 

27.2 

2.3 

5 

15 

36 

108 

9 

59.72 

5.62 

3.2 

6 

9 

51 

179 

10 

69.2 

34.7 

2.6 

5 

15 

12 

57 

Average  over  10  events 

16 

26 

68 

Table  2 
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Figure  Captions 


Fig.  1.  NORESS  noise  correlations  versus  interstation  separation  for  the  two  frequency 
bands  1-2  Hz  and  3-4  Hz.  The  noise  segment  used  is  30  s  long  and  taken  at  05:l'5  GMT  on 
day  323  of  1985.  Mean  values  and  standard  deviations  within  100  m  distance  intervals  are 
plotted  on  top  of  the  population,  except  for  short  and  long  distances,  where  the  number 
of  correlation  values  is  low. 

Fig.  2.  The  geometry  of  the  NORESS  array.  The  instrument  at  the  center  is  denoted 
AO.  The  other  sensors  are  arranged  in  four  concentric  rings,  the  A-ring,  B-ring,  C-ring 
and  D-ring.  The  type  of  instrumentation  at  each  of  the  sensor  sites  is  given  in  the  legend. 

Fig.  3.  NORESS  transfer  functions  (velocity  sensitivity)  for  (from  left  to  right)  the  long 
period,  intermediate  period,  short  period  and  high  frequency  band. 

Fig.  4.  The  figure  shows  data  from  the  vertical  short  period  sensor  at  NORESS  site  AO 
for  an  explosion  in  the  White  Sea  region  of  the  USSR  on  18  July  1985.  (ISC  solution: 
origin  time  21.14.57.7,  epicenter  at  65.96°,  40.86°E,  depth  0  km  and  mj,  5.1).  Also  shown 
are  broad  band  frequency- wavenumber  spectra  computed  for  the  short  data  windows  (in¬ 
dicated  on  top  of  the  trace)  centered  around  the  onsets  of  the  Pn,  Sn  and  Lg  arrivals. 

Fig.  5.  NORESS  data  for  a  small  event  in  western  Norway  at  epi central  distance  350  km. 
The  various  regional  phases  are  enhanced  by  frequency  filtering  in  different  passbands,  as 
explained  in  the  text.  AOZ  denotes  the  vertical  component  of  the  short  period  sensor  at 
the  array  center,  whereas  AOE  is  the  horizontal  east-west  component  at  the  same  site. 

Fig.  6.  NORESS  event  processor  plot  for  a  mining  explosion  in  Estonia.  The  plot  shows 
data  for  the  vertical  channels  of  AO  and  all  instruments  in  the  B-  and  C-rings.  The  two 
bottom  traces  are  tlie  beams  on  which  the  Pn  (beam  N073)  and  Lg  phase  (beam  NV04) 
were  detected.  Note  that  the  incoherent  beam  is  plotted  here  as  a  coherent  beam,  but  with 
the  steering  delays,  filters  and  configuration  as  defined  for  the  incoherent  beam.  Vertical 
bars  across  tlie  panel  indicate  detected  phases  lh.at  are  used  in  the  TTAZLOC  scheme 
(Bratt  and  Bache,  1988)  to  locate  the  event.  Bulletin  information  is  given  in  the  upper 
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line  on  top  of  the  traces.  The  second  line  contains  information  pertinent  to  the  Pn  phase. 

Fig.  7.  The  figure  shows  168  NORESS  noise  suppression  spectra  taken  hourly  during  a 
one-week  period  in  July  1986.  The  subgeometry  used  comprises  the  AO,  C-  and  D-ring 
sensors.  The  horizontal  line  indicates  the  i/iV  (N  =  17  in  this  case)  suppression  level. 

Fig.  8.  P-phase  detection  statistics  for  ARCESS  for  regional  events  in  the  distance  range 
800-1200  km,  using  the  Helsinki  bulletin  as  a  reference.  The  upper  part  of  the  figure  shows 
the  distribution  of  events  by  magnitude  with  detected  events  corresponding  to  the  hatched 
columns.  The  bottom  part  of  the  figure  shows  the  estimated  detection  probability  curve 
as  a  function  of  magnitude,  with  the  observed  detection  percentages  marked  as  asterisks. 
The  stippled  curves  mark  the  90  per  cent  confidence  limits. 

Fig,  9.  NORESS  and  ARCESS  data  for  the  sensor  at  Cl  for  a  nuclear  test  at  Novaya 
Zemlya  on  December  4,  1988.  Data  are  shown  both  unfiltered  (upper  trace)  and  filtered 
in  the  band  5-10  Hz  (lower  trace)  for  both  arrays. 

Fig.  10,  The  figure  shows  the  results  of  the  narrow  band  (left)  and  broad  band  (right) 
f-k  estimator  applied  to  the  Pn,  Sn  and  Lg  phases  observed  at  NORESS  from  a  suite  of 
10  chemical  explosions  at  a  dam  construction  site  in  southern  Norway.  The  processing 
frequencies  for  the  narrow  band  method  were  7  Hz  for  Pn  and  4  Hz  for  Sn  and  Lg.  For 
the  broad  band  method  estimates  are  given  for  the  frequency  intervals  5.25-8.75  Hz  for 
Pn  and  3-5  Hz  for  Sn  and  Lg. 

Fig.  11.  iilustralioii  of  variation  of  relative  importance  of  tiie  plia.se.s  Sn  and  Lg  for  .six 
events  with  locations  as  indicated  in  the  map.  The  standard  group  velocities  of  4.5  and 
3.5  km/s,  commonly  assigned  to  Sn  and  Lg,  respectively,  are  marked  by  dashed  lines. 
The  upper  three  traces  cover  the  distance  interval  480-550  km,  while  the  lower  three 
traces  correspond  to  epicentral  distances  in  the  range  1225-1320  km.  The  location  of 
the  NORSAR  array  is  denoted  by  a  ring  on  the  map,  and  the  traces  are  from  NORSAR 
seismometer  02B01.  The  data  are  bandpass  filtered  1-5  Hz.  The  reduction  velocity  is  8.0 
km/s. 
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Fig.  12.  Azimuth  residuals  (AR  =  estimated  azimuth  minus  ‘true’  azimuth,  in  degrees) 
for  Pn  phases  detected  at  NORESS  during  1985-88,  for  events  rei)orted  in  the  Helsinki 
bulletin.  Sec  the  text  for  further  explanation. 

Fig.  13.  The  map  shows  the  location  of  the  three  regional  arrays  NORESS  and  ARCESS 
in  Norway  and  FINESA  in  Finland,  as  v'ell  as  the  location  of  the  ten  events  used  in  the 
TTAZLOC  location  estimation  experiment. 
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ARCESS  P -PHASE  DETECTION 

REFERENCE  HELSINKI  BULLETIN  JAN-MAR  88 

T.ENINGRAD  REGION  ,  DISTANCE  800-1200  KM 


Fig.  10 
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Dr.  Stephen  Bratt 
Center  for  Seismic  Studies 
1300  North  17th  Street 
Suite  1450 

Arlington,  VA  22209 

Michael  Browne 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Mr.  Roy  Burger 
1221  Serry  Road 
Schenectady,  NY  12309 


Dr.  Robert  Burridge 
Schlumberger-Doll  Research  Center 
Old  Quarry  Road 
Ridgefield,  CT  06877 


Dr.  W.  Winston  Chan 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314-1581 


Dr.  Theodore  Cherry 
Science  Horizons,  Inc. 

710  Encinitas  Blvd.,  Suite  200 
Encinitas,  CA  92024  (2  copies) 


Prof.  Jon  h.  Claerbout 
Department  of  Geophysics 
Stanford  University 
Stanford.  CA  94305 


Prof.  Robert  W.  Clayton 

Seismological  Laboratory 

Division  of  Geological  &  Planetary  Sciences 

California  Institute  of  Technology 

Pasadena,  CA  91125 

Prof.  F.  A.  Dahlen 

Geological  and  Geophysical  Sciences 

Princeton  University 

Princeton,  NJ  08544-0636 


5 


Prof.  Adam  Dziewonski 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  St 
Cambridge,  MA  02138 

Prof.  John  Ebel 

Department  of  Geology  &  Geophysics 
Boston  College 
Che.stnut  Hill,  MA  02167 


Eric  Fielding 
SNEE  Htdl 
INSTOC 

Cornell  U  in’w.  sity 
Ithaca,  N'i  14853 

Prof.  Donald  Forsyth 
Department  of  Geological  Sciences 
Brown  University 
Providence,  RI  02912 


Dr.  Cliff  Frolich 
Institute  of  Geophysics 
8701  North  Mopac 
Austin,  TX  78759 


Dr.  Anthony  Gangi 
Texas  A&M  University 
Department  of  Geophysics 
College  Station,  TX  77843 


Dr.  Freeman  Gilbert 
IGPP,  A-025 
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P.O.  Box  51 

N-2007  KjeUer,  NORWAY 


Dr.  Michel  Campillo 
Observatoirc  de  Grenoble 
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